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依据 单 胺 缺陷 假说 (monoamine deficiency hypothesis) ， 多 巴 胺 系统 功能 失调 是 引发 抑 
郁 的 重要 原因 (Belmaker & Agam, 2008) 。 多 巴 胺 系统 功能 受到 多 巴 胺 代谢 、 转 运 和 传导 多 
个 环节 的 共同 影响 (Opmeer, Kortekaas, & Aleman, 2010) 。 此 外 ， 越 来 越 多 的 研究 显示 中 脑 
一 边缘 多 巴 胺 通路 (mesolimbic pathway) 和 中 脑 一 皮质 多 巴 胺 通路 〈mesocortical pathway) 
功能 缺陷 与 抑 邦 密切 相关 (Dunlop & Nemeroff, 2007) 。 基 于 此 , 在 众多 多 巴 胺 系统 基因 中 ， 
既 参与 调节 多 巴 胺 代谢 、 转 运 和 传导 ， 又 表达 于 中 脑 一 边缘 、 中 脑 一 皮质 多 巴 胺 通路 中 的 基 
因 位 点 受到 了 研究 者 的 广泛 关注 。 壁 如 ， 儿 茶 酚 胺 氧 位 甲 基 转 移 酶 基因 (catechol-0- 
methyltransferase, COMT) 、 多 巴 胺 转运 体 基 因 (dopamine transporter, D471) 及 D2 型 多 巴 


胺 受 体 基 因 (dopamine receptor D2, DRD2) (Cao etal., 2018; Lin et al., 2017; Pinsonneault et 


al., 2011) 。 
COMT 基因 编码 的 COMT MERER BEARS A A LRR AERA ZEREA AS, 
_ 降低 突 触 间隙 中 的 多 巴 胺 浓度 进而 导致 抑郁 。 在 众多 COMT 基因 多 态 性 中 ， 多 数 研究 者 关 
= 注 rs4680 (Vall 58Met) 多 态 性 与 抑郁 等 心理 病理 问题 的 关联 , 但 是 研究 显示 在 亚洲 样本 中 ， 
rs6267 多 态 性 与 COMT 酶 活性 的 关联 强 于 rs4680 多 态 性 (Lee et al.,2005) 。 基 于 此 ， 本 研 
究 重点 关注 位 于 22 号 染色 体 q11.2 区 的 rs6267 多 态 性 ,在 该 多 态 性 中 ,与 T 等 位 基因 相 比 ， 
G 等 位 基因 与 更 高 COMT 酶 活性 有 关 (Leeetal.,2005) ， 即 G 等 位 基因 具有 较 高 的 多 巴 胺 
= 代谢 率 进而 会 导致 较 低 的 多 巴 胺 浓度 。Lin (2017) 的 研究 显示 ，rs6267 多 态 性 与 帕 金 森 
一 病人 的 抑郁 症状 显著 相关 。 
x DATI 基因 编码 的 多 巴 胺 转运 体能 够 将 突 触 间隙 的 多 巴 胺 转运 再 摄取 到 突 触 前 膜 ， 调 节 
脑 内 多 巴 胺 含量 进而 影响 抑郁 。 现 有 研究 主要 考察 了 DATI VNTR 多 态 性 ， 但 是 在 何 种 等 位 
= 基因 COR BK 10R) 与 较 高 的 DAT 表达 活性 有 关上 仍 存 在 分 歧 (e.g., Costa, Riedel, Müller, 
一 Miller, & Ettinger, 2011; Heinzetal.,2000) 。 相 比 之 下 ， 位 于 5 号 染色 体 p15.33 区 的 rs27072 
多 态 性 表达 活性 更 为 稳定 , 在 体内 和 体外 实验 中 一 致 发 现 T 等 位 基因 比 CC 基因 型 具有 更 高 
的 表达 水 平 (Pinsonneaultetal., 2011) 。 有 关 双 相 障 碍 的 研究 显示 , 在 多 个 重复 验证 样本 中 ， 
相 比 其 他 基因 位 点 ，rs27072 多 态 性 与 双 相 障碍 存在 稳定 关联 (Pinsonneault et al., 2011) 。 
DRD2 基因 编码 的 D2 型 多 巴 胺 受 体 通过 影响 中 脑 边缘 系统 的 多 巴 胺 能 活性 从 而 影响 众 
患 抑郁 的 风险 。rs1799978 多 态 性 位 于 11 号 染色 体 q2.2-2.3 区 ， 是 DRD2 基因 启动 子 区 的 一 
种 功能 性 多 态 性 位 点 ,在 -241 位 置 处 腺 味 叭 (A) 置 换 鸟 味 吟 (G) 形 成 两 种 等 位 基因 (Doehring， 
Kirchhof, & Létsch, 2009) ， 影 响 D2 受 体 的 表达 (Zhang & Malhotra, 2011) 。 尽 管 尚未 有 关 
于 rs1799978 等 位 基因 功能 的 直接 证 据 , 但 是 相关 研究 提示 A 等 位 基因 可 能 与 更 低 的 多 巴 胺 
受 体 活性 (Doehring etal., 2009) 及 更 高 的 抑郁 风险 有 关 《〈Cao et al,2018) 。 有 具体 而 言 ， 有 
关 利 培 酮 (通过 与 D2 受 体 结合 产生 拷 抗 作用 以 减少 抑郁 、 精 神 分 裂 症 状 ) 的 药物 研究 显示 ， 
相 比 G 等 位 基因 ，A 等 位 基因 或 AA 基因 型 携带 者 在 治疗 后 表现 出 更 好 的 疗效 (Doehring et 
al.,2009) 。 这 提示 ，A 等 位 基因 可 能 与 较 低 的 D2 受 体 数量 相关 ， 较 低 的 D2 受 体 数 量 可 以 
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使 得 利 培 酮 与 受 体 结合 率 更 高 ， 从 而 起 到 更 好 的 疗效 。 此 外 ， 相 关 研究 显示 在 亚 裔 样本 中 ， 
相 比 其 他 DRD2 多 态 性 ，rs1799978 多 态 性 与 更 好 的 利 培 酮 疗效 有 关 (Xing et al., 2007) ， 
并 且 与 抑郁 的 关联 更 强 (Cao etal., 2018) 。 
尽管 单 基因 研究 为 理解 抑郁 的 发 生机 制 与 个 体 差异 提供 了 重要 启示 , 但 是 单 基因 研究 的 
解释 率 往往 极 低 。 造 成 微弱 的 遗传 解释 率 的 原因 之 一 可 能 是 研究 忽视 了 抑郁 的 多 基因 遗传 结 
构 (Belsky & Pluess, 2009) 。 尤 其 是 ， 研究 显 示 COMT ÆA, DATI 基因 和 DRD2 基因 可 能 
并 非 独立 作用 于 抑郁 , 而 是 存在 联合 效应 ,如 早期 的 动物 研究 发 现 , 在 DATI FEAR EE 
COMT iit St LEAR BR Bis 400%, D2 受 体 表达 水 平 下 降 (Giros, Jaber, Jones, Wightman, & 
Caron, 1996) 。 这 表明 COMT 基因 、D477 基因 和 DRD2 基因 间 可 能 存在 联合 效应 ， 共 同调 
节 多 巴 胺 系统 的 功能 。 
近期 , 研究 者 通过 累加 个 体 携带 的 易 感 等 位 基因 条 数 测量 多 基因 联合 效应 , 即 多 基因 累 
— 加 得 分 (multilocus genetic profile score, MGPS) 。 多 基因 累加 得 分 方法 为 抑郁 的 遗传 机 种 
究 提供 了 一 种 新 的 框架 ， 人 迄今 已 有 研究 显示 多 基因 效应 比 任何 单一 位 点 的 效应 更 强 
(Pearson-Fuhrhop et al., 2014; Vrshek-Schallhor et al., 2015) 。 璧 如 ，Pearson-Fuhrhop 等 人 
《2014) 依 据 单 胺 缺陷 假说 ,考察 了 影响 多 巴 胺 代谢 、 转 运 和 传导 功能 的 5 种 DA 基因 (COMT、 
DATI, DRDI, DRD2 和 DRD3) 的 多 基因 累加 得 分 与 抑郁 症 的 关联 ， 结 果 显 示 无 论 是 在 抑 
郁 患 者 组 还 是 健康 对 照 组 ， 个 体 携带 的 风险 等 位 基因 ( 低 活性 多 巴 胺 等 位 基因 〉 数量 越 多 ， 
其 抑郁 水 平 越 高 ， 且 相 比 单个 基因 位 点 ,多 基因 累加 得 分 的 解释 率 更 高 。 这 提示 多 基因 累加 
方法 能 够 提供 更 丰富 的 遗传 信息 并 且 增 加 遗传 解释 力 。 此 外 ,在 统计 分 析 方 面 ， 多 基因 累加 
得 分 可 作为 一 种 连续 变量 , 在 检验 基因 x 环境 交互 作用 时 更 具有 统计 优势 (Dunn et al., 2011). 
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12 多 基因 累加 得 分 与 环境 因素 的 交互 作用 

理论 和 实证 研究 显示 , 多 巴 胺 系统 多 基因 不 仅 联 合影 响 个 体 乌 患 抑 郁 的 风险 , 并 且 可 能 
调节 个 体 对 环境 的 敏感 性 。 壁 如 ,生物 社会 发 展 模型 (biosocial developmental model, BDM) 
指出 , 多 巴 胺 舍 量 下 降 会 引发 青少年 行为 趋向 系统 功能 缺陷 , 增强 个 体 在 面 对 不 利 环 境 时 的 
敏感 性 ， 导 致 个 体 乌 患 内 外 化 问题 的 风险 提高 (Beauchaine, Gatzke-Kopp & Mead, 2007) 。 
前 额 叶 晚熟 理论 认为 , 边缘 系统 一 前 额 叶 连通 性 影响 个 体 对 环境 的 敏感 性 , 进而 与 抑郁 情绪 
密切 相关 (Andersen & Teicher, 2008) 。COMT、DAT1 和 DRD2 基因 的 主要 表达 脑 区 分 别 集 
中 于 大 脑 前 额 叶 区 域 (Matsumoto et al., 2003)、 中 脑 纹 状 体 及 边缘 系统 区 域 (Lewis, Melchitzky, 
Sesack, Whitehead, & Sampson, 2001; Noble, Gottschalk, Fallon, Ritchie, & Wu, 1997) 。 因 此 ， 
作用 于 边缘 系统 一 前 额 叶 通路 的 三 种 基因 可 能 联合 调节 个 体 对 环境 的 敏感 性 进而 影响 抑 邦 
的 发 生发 展 。 实 证 研究 也 为 多 巴 胺 系统 多 基因 累加 得 分 与 环境 的 交互 作用 提供 了 初步 的 证 
据 。 如 Davies, Cicchetti 和 Hentges (2015) 的 研究 显示 ， 在 经 历 母 亲 消 极 教养 行为 时 ， 多 
基因 累加 得 分 较 高 的 青少年 比 多 基因 累加 得 分 较 低 的 青少年 问题 行为 更 多 。 关 于 上 腹 侧 纹 状 体 
激活 水 平 的 研究 显示 ， 在 面 对 奖 赏 信息 反馈 时 ， 个 体 携带 的 低 活性 多 巴 胺 系统 等 位 基因 
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(DRD4 Æ, DATI 基因 、DRD2 基因 和 COMT 基因 ) 数量 越 多 ， 则 表现 出 更 低 的 纹 状 体 
反应 性 (Nikolova, Ferrell, Manuck, & Hariri, 2011) 。 

然而 ， 迄 今 少数 几 项 多 巴 胺 系统 多 基因 累加 得 分 x 环 境 研究 主要 关注 不 利 的 家 庭 环 境 因 
素 , 如 消极 教养 (Davies etal., 2015)、 产 前 不 利 环境 (Bischoffet al., 2017) . HUE FE (Coley, 
Sims, & Carrano, 2017) 等 ， 尚 未 有 研究 关注 积极 环境 的 作用 。 在 众多 的 家 庭 环 境 因 素 中 ， 母 
亲 积 极 和 消极 教养 行为 是 影响 青少年 抑郁 的 重要 预测 源 。 研 究 显 示 ， 母 亲 的 积极 教养 行为 ， 
如 温情 、 支持 性 的 教养 行为 会 促进 青少年 积极 的 自我 知觉 和 较 好 的 社会 情绪 功能 , 进而 降低 
抑郁 的 风险 (Olino et al., 2016; Wang et al., 2016) 。 相 反 ， 消 极 教养 行为 ， 如 拒绝 、 严 厉 等 
会 导致 青少年 形成 较 差 的 自我 意识 和 自尊 ， 继 而 导致 较 高 的 抑郁 风险 (Wang et al., 2016) 。 
因此 ， 本 研究 的 目的 之 一 是 综合 运用 积极 和 消极 教养 指标 ， 考 察 COMT 基因 、D477 基因 和 
DRD2 基因 多 基因 累加 得 分 与 环境 的 潜在 交互 作用 。 


一 1.3 “素质 -压力 ”模型 与 “不 同 易 感性 ”模型 
‘ 多 基因 累加 得 分 与 环境 的 交互 作用 模式 可 能 被 两 种 竞争 假设 所 解释 : “素质 -压力 ” 模 
型 与 “不 同 易 感性 ”模型 。 前 者 假设 一 些 心理 疾患 或 行为 问题 的 产生 是 由 于 个 人 携带 的 风险 
基因 和 不 利 环境 的 综合 结果 (Monroe & Simons, 1991) ; 而 后 者 认为 ， 所 谓 的 风险 基因 事实 
上 是 一 种 遗传 可 塑性 ， 更 具 可 塑性 的 个 体 对 消极 和 积极 环境 均 更 加 敏感 (Belsky & Pluess, 
2009) 。 基 于 这 两 种 竞争 模型 ， 研 究 者 对 多 基因 累加 得 分 本 质 上 是 一 种 “累加 风险 得 分 ”还 
是 “累加 可 塑性 得 分 ” 亦 存 在 纷争 。 

然而 问题 是 , 当前 的 多 巴 胺 系统 多 基因 累加 得 分 研究 中 所 测量 的 环境 指标 主要 局 限于 消 
极 环境 ， 尤其 是 当 其 取 值 范围 有 限时 , 采用 传统 回归 方法 可 能 无 法 反映 交互 作用 的 全 貌 ， 其 
所 验证 的 “素质 -压力 ”模型 可 能 反映 了 “不 同 易 感性 ”模型 的 一 部 分 (Roisman etal., 2012)。 
基于 此 ， 研 究 者 提出 了 更 为 准确 的 、 新 兴 的 统计 检验 方法 一 一 再 参数 化 回归 模型 Cre- 


oO parameterized regression model) ， 能 够 从 理论 假设 出 发 检验 “基因 x 环 境 ” 交 互 作用 的 形式 
(Widaman et al., 2012) 。 因 此 ， 本 研究 的 第 二 个 目的 是 采用 新 兴 的 再 参数 化 回归 模型 验证 
两 种 竞争 假设 。 


综 上 ， 本 研究 采用 多 基因 研究 视角 ， 考 察 多 巴 胺 系统 基因 累加 得 分 (COMT、DATI 和 
DRD2) 与 母亲 教养 行为 对 青少年 抑郁 的 影响 ， 主 要 探讨 以 下 问题 : (1) 多 基因 累加 得 分 和 
母亲 教养 行为 对 青少年 抑郁 的 交互 作用 ; D 对 多 基因 累加 得 分 与 母亲 教养 行为 影响 抑郁 
的 交互 作用 形式 (“素质 -压力 ”模型 vs.“ 不 同 易 感 性 ”模型 ) 进行 检验 。 


2 WAI 
2.1 被 试 

被 试 来 自 国 内 一 项 大 型 追踪 项 目 。 该 项 目 最 初 在 山东 省 济南 市 的 14 所 小 学 整 群 抽取 三 
年 级 小 学 生 2261 名 ， 在 每 年 第 二 学 期 进行 一 次 追踪 测 查 。 研 究 显 示 ， 在 童年 期 ， 抑 郁 的 发 
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生 率 和 严重 程度 较 低 (Cole etal., 2002) ， 但 是 六 、 七 年 级 〈 约 12713 岁 ) 是 抑郁 水 平 迅速 
上 升 的 关键 期 ， 而 之 后 抑郁 水 平 的 增长 速度 减缓 并 趋 近 于 成 人 水 平 〈Cole et al., 2002; Ferro, 
Gorter, & Boyle, 2015; Natsuaki, Biehl, & Ge, 2009) ， 因 此 本 研究 选择 在 青少年 抑郁 快速 上 升 
的 年 龄 阶段 进行 研究 。 鉴 于 工作 量 和 研究 经 费 问题 ， 以 各 学 校 班级 为 单位 随机 抽取 部 分 被 试 

(1086 A) 进行 遗传 数据 收集 。 根 据 Duncan fil Keller (2011) 提出 的 遗传 研究 统计 检验 力 
与 样本 量 需 求 ， 当 估计 遗传 效应 量 在 1% 左 右 时 ， 为 获得 80% 以 上 的 统计 检验 力 ， 样 本 量 要 
保证 在 600 人 以 上 , 本 研究 被 试 量 满足 最 低 样本 量 要 求 。 具有 遗传 数据 的 被 试 与 原始 样本 在 
基本 研究 变量 性别、 年 龄 、 抑 郁 、 母 亲 教 养 等 ) 上 均 无 显著 差异 (性 别 : x?=2.19, df= 1,p 
= 0.14; 母亲 受 教育 水 平 : X= 2.83, df= 2,p= 0.24; 父亲 受 教育 水 平 :X= 2.61, 4f= 2,p= 0.27; 
家 庭 月 收入 : 72 = 4.87, df = 2, p = 0.09; 年 龄 : te18s) = 0.12, p = 0.91; T1 抑郁 : teri) = 1.53, p = 
0.13; T2 抑郁 : teon = 0.80, p = 0.43; 积极 教养 : t3206) =—0.91, p = 0.36; 消极 教养 : t3206) = 0.76, 
p=0.45) 。 同 时 ， 为 避免 人 口 分 层 偏 差 ， 本 研究 仅 在 汉族 被 试 (1052 人 〉 中 进行 分 析 ， 汉 
族 被 斌 与 少数 民族 被 试 在 基本 研究 变量 上 亦 无 显著 差异 (性 别 : x? = 0.00, df= 1,p = 0.98, 年 
龄 : taos3) = 1.47, p = 0.14; COMT: X=3.71, df= 2, p = 0.168; DAT1: x? = 0.55, df= 2, p = 0.76; 


ant 
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DRD2: 72 = 1.19, df= 2, p = 0.55; T1 抑 邦 : tao72) = 0.06, p = 0.96; T2 抑郁 : taom = 0.49, p = 0.63; 
积极 教养 : tq073)= 0.32,p=0.75; 消极 教养 :xiao7za)= 0.03, p=0.97) 。 最 终 ， 本 研究 被 试 共 1052 
人 ， 六 年 级 时 平均 年 龄 为 12.31 + 0.37 岁 ， 男 生 524 人 (49.8%) ， 女 生 528 人 (50.2%) 。 
母亲 受 教育 水 平 在 初中 及 以 下 者 占 11.7%， 高 中 或 中 专 者 占 26.5%， 大 专 以 上 者 占 61.8%。 
父亲 受 教育 水 平 在 初中 及 以 下 者 占 8.3%， 高 中 或 中 专 者 占 20.6%， 大 专 以 上 者 占 71.0%。 家 
庭 月 收入 低 于 3000 元 的 占 23.2%，3000 ~ 6000 元 的 占 47.9%， 高 于 6000 元 的 占 28.8%. 


22 ”研究 工具 
2.2.1 ”母亲 教养 行为 

采用 中 文 版 儿童 养育 实践 问卷 (Child-Rearing Practices Report) (Chen, Bian, Xin, Wang, 
& Silbereisen, 2010) 测量 母亲 教养 行为 ， 该 问卷 由 母亲 报告 。 采 用 温情 〈4 个 项 目 ， 如 “我 
以 温和 、 亲 切 的 态度 和 孩子 说 话 ”) 、 引 导 (4 个 项 目 ， 如 “ 当 我 的 孩子 遇 到 什么 问题 时 ， 
我 鼓励 他 〈 她 ) 把 它 说 出 来 ”) 2 个 维度 作为 积极 教养 行为 的 指标 ， 采 用 拒绝 (4 个 项 目 ， 
如 “如 果 我 的 孩子 不 来 烦 我 ， 我 就 不 会 理 皮 她 ”) 作为 消极 教养 行为 指标 。 问 卷 采用 5 点 记 
分 (0 表示 完全 不 符合 ，4 表示 完全 符合 ) ， 平 均 得 分 越 高 ， 则 该 种 教养 行为 越 多 。 本 研究 
中 积极 和 消极 教养 行为 的 Cronbach’s a 系数 分 别 为 0.83 和 0.54。 验 证 性 因素 分 析 显 示 二 因 
子 模型 拟 合 良好 (x? = 223.58, df= 53, RMSEA = 0.06, CFI = 0.92, TLI = 0.90) 。 


2.2.2 ”青少年 抑郁 
采用 儿童 抑郁 量 表 (Children’s Depression Inventory, CDI; Kovacs, 1992) 测量 青少年 抑 
郁 症状 。 该 量 表 共 包 含 27 个 项 目 ， 要 求 青少年 报告 近 两 周 内 的 抑郁 症状 ， 分 别 采 用 0、1、 


2 记分 ， 计 算 青 少年 的 抑郁 均 分 ， 得 分 越 高 则 表明 青少年 的 抑郁 症状 越 多 。 该 量 表 在 正常 青 
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少年 群体 中 应 用 广泛 (Zhang et al., 2015) ， 具 有 良好 的 测量 学 特性 。 本 研究 中 2 个 时 间 点 
青少年 抑郁 的 Cronbach’s a 系数 均 为 0.88。 验 证 性 因素 分 析 显 示 T1 抑郁 G2 = 699.34, df= 
324, RMSEA = 0.03, CFI = 0.95, TLI = 0.95) 和 T2 抑郁 (X = 736.72, df= 324, RMSEA = 0.04, 
CFI = 0.96, TLI = 0.96) 的 单 因子 模型 均 拟 合 良 好 。 


2.3 ”研究 程序 

研究 经 过 本 单位 伦理 委员 会 审核 批准 。 首 先 , 项 目 组 将 问卷 测 查 、 唾 液 样品 采集 、DNA 
提取 与 分 型 流程 等 信息 告知 施 测 学 校 、 青 少年 父母 及 青少年 本 人 , 获得 三 方 的 知情 同意 后 进 
行 了 数据 收集 。 其 次 ， 所 有 施 测 程序 均 由 经 过 严格 培训 的 研究 生 主 试 完成 。 
2.3.1 问卷 施 测 流程 

首先 ， 以 班级 为 单位 , 采用 现场 问卷 调查 的 方式 测评 青少年 抑郁 症状 , 每 个 施 测 班 级 由 
2 名 经 过 严格 培训 并 具有 丰富 施 测 经 验 的 研究 生 担 任 主 试 , 施 测 结束 后 问卷 当场 收回 。 其 次 ， 
现场 测评 后 ,由 青少年 将 装 有 母亲 教养 问卷 的 信封 带 回 家 交 由 母亲 填写 , 第 二 天 带 回 交 给 班 
主任 老师 统一 回收 。 
2.3.2 ”唾液 样本 采集 、DNA 提取 与 分 型 流程 

首先 ， 以 班级 为 单位 采集 青少年 的 唾液 样本 ， 保 证 每 人 采集 2ml 以 上 。 采 集 唾液 样品 
后 ， 利 用 Sequenom (San Diego, CA, USA) 芯片 基质 辅助 激光 解吸 /电离 飞行 时 间 (MALDI- 
TOF) 质谱 平台 对 rs6267、rs27072、rs1799978 多 态 性 位 点 进行 DNA 提取 与 基因 分 型 。 PCR 


引物 分 别 为 ，COMT 基因 : forward 5’°-ACGTTGGATGTAGGTGTCAATGGCCTCCAG-3’, 


reverse 5’°-ACGTTGGATGTCATGGGTGACACCAAGGAG-3’; DATI 基因 : forward 5’- 
AGAACACAGTGCCCCTGGG-3’, reverse 5’°-AAAAACGTCTAACTTCATGCTGTCTG-3’; DRD2 


基 因 : forward 5’-AGGAGCTGGAGATGGAGATGCT-3’ i reverse 5- 


ATGCCCATTCTTCTCTGGTTTGGC-3’. PCR 反应 条 件 : 94°C 15 min; 94°C 20 s, 56°C 30 s, 
72C 1 min, tt 45 个 循环 ， 最 终 72°C 3 min。 随 后 经 历 单 碱 基 延伸 反应 ， 基 因 分 型 采用 
MassARRAY Typer 3.4 软件 系统 进行 分 析 。 本 研究 所 使 用 的 检测 平台 和 技术 具有 较 高 的 可 靠 
性 (基因 分 型 有 效率 > 97%) 。 


2.4 数据 处 理 与 分 析 

首先 ， 对 3 个 基因 多 态 性 进行 Hardy-Weinberg 平衡 检验 并 考察 不 同 基因 间 的 关联 。 其 
次 ， 对 3 个 基因 的 线性 基因 效应 ! 和 等 基因 效应 ?进行 检验 ， 确 定 线 性 基因 编码 方式 是 否 合适 
以 及 单个 基因 是 否 存在 主导 效应 (Stocker et al.,2017) 。 该 部 分 统计 分 析 共 检验 三 个 回归 模 


1 线性 基因 效应 : 基因 编码 方式 影响 遗传 研究 结果 (Aliev, Latendresse, Bacanu, Neale, & Dick, 2014) 。 迄 今 ， 遗 传 研究 可 能 采 
显 性 编码 (0= 不 携带 风险 等 位 基因 ; 1= 携 带 1 或 2 个 风险 等 位 基因 〉、 隐 性 编码 (0= 携 带 0 或 1 个 风险 等 位 基因 ; 1= 携 带 2 
个 风险 等 位 基因 ) 和 线性 编码 (0= 不 携带 风险 等 位 基因 ; 1= 携 带 1 个 风险 等 位 基因 ; 2= 携 带 2 个 风险 等 位 基因 ) 。 多 基因 累加 
风险 研究 通常 采用 线性 编码 方式 而 忽视 遗传 基因 在 功能 上 是 否 符合 线性 模式 ， 因 此 有 必要 对 不 同 基因 型 的 线性 效应 进行 检验 。 
2 等 基因 效应 : 多 基因 累加 风险 研究 通常 将 多 个 基因 的 效应 进行 累加 ， 但 是 无 法 排除 单个 效应 极 大 的 基因 位 点 主导 多 基因 效应 
情况 ， 此 外 若 忽 视 不 同 基 因 的 功能 ， 随 机 的 基因 编码 方式 可 能 存在 相似 基因 间 的 效应 编码 相反 ， 引 发 多 基因 间 的 效应 相互 抵 
因此 ， 有 必要 对 多 基因 的 编码 方式 进行 效应 大 小 和 方向 的 等 同性 检验 。 
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型 : 分 解 模型 〈disaggregated model) 包括 6 个 基因 主 效应 (3 个 基因 的 显 性 和 隐 性 基因 编 
码 ) 、6 个 基因 x 环 境 交 互 效应 〈 即 3 个 基因 的 显 性 和 隐 性 基因 编码 与 教养 的 交互 作用 ) 、 
教养 和 控制 变量 的 主 效应 《性别 和 早期 抑郁 ) 共 15 个 参数 估计 ; 线性 基因 效应 模型 〈linear 
gene effect model) 用 于 检验 线性 基因 编码 是 否 合适 , 将 同一 个 基因 的 显 性 和 隐 性 编码 回归 权 
EREE (6 个 限定 ， 包 括 基因 主 效应 和 基因 x 环境 交互 效应 ， 人 允许 不 同 基因 间 的 效应 存 
在 差异 ) ; 等 基因 效应 模型 (equal gene effect) 将 3 个 基因 的 效应 限定 相等 (4 MRE, E 
括 基因 主 效应 和 基因 x 环境 交互 效应 ) 。 线 性 基因 效应 和 等 基因 效应 模型 的 优势 在 于 可 以 获 
得 模型 限定 的 估计 : 若 线性 基因 效应 相 比分 解 模 型 有 显著 的 解释 率 变化 则 表明 至 少 一 个 基因 
存在 非 线 性 效应 ， 则 该 基因 不 适合 0、1、2 的 线性 编码 方式 ， 需 要 调整 基因 编码 方式 ; 若 等 
基因 效应 模型 比 线性 基因 效应 模型 有 明显 的 解释 率 变化 则 表明 3 个 基因 的 预测 强度 存在 明 
显 差异 , 至 少 一 个 基因 存在 不 同 的 预测 强度 或 预测 方向 , 可 能 存在 单个 基因 的 主导 效应 或 基 
因 间 效应 相互 抵消 ， 不 适合 进行 累加 基因 得 分 研究 。 第 三 , 采用 相关 分 析 考 察 主要 研究 变量 
间 的 关系 。 第 四 , 采用 分 层 线性 回归 考察 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 抑郁 的 影 
响 ， 为 避免 多 重 共 线 性 ， 多 基因 累加 得 分 和 母亲 教养 行为 均 进行 标准 化 处 理 ， 若 基因 x 环 境 
交互 作用 显著 则 进一步 进行 简单 斜率 检验 。 同 时 ， 将 性 别 和 早期 抑郁 症状 作为 控制 变量 ， 控 
Hl) T1 抑郁 本 质 上 考察 了 基因 与 教养 对 Tl 到 T2 抑郁 变化 的 影响 。 第 五 ， 采 用 再 参数 化 回归 
(Widaman et al., 2012) 检验 交互 作用 符合 何 种 理论 模型 (“素质 -压力 ”和 “不 同 易 感性 ”) 。 
最 后 ， 为 验证 结果 的 可 靠 性 进行 敏感 性 分 析 。 敏 感性 分 析 分 五 个 方面 进行 ，(1) 将 每 个 基 
因 与 母亲 教养 行为 对 抑郁 的 交互 作用 进行 检测 , 比较 单 基因 和 多 基因 结果 的 显著 性 差异 , (2) 
将 多 基因 累加 得 分 划分 为 低 、 中 、 高 三 组 ， 进 行 分 组 回归 分 析 验 证 多 基因 x 环境 交互 结果 。 
i (3) 对 基因 x 基因 交互 效应 进行 检验 , 排除 基因 间 交 互 效应 对 多 基因 累加 效应 的 干扰 。(4) 
= 参照 既 有 研究 文献 (Cao et al., 2018) ， 将 消极 教养 的 项 目 全 部 进行 反 向 计 分 ， 随 后 将 积极 
PT 和 消极 教养 ( 反 向 计 分 后 ) 的 项 目标 准 化 后 计算 均值 作为 母亲 教养 敏感 性 得 分 , 重复 上 述 分 
析 步 骤 ， 用 以 验证 结果 的 可 靠 性 。 为 控制 1 型 错误 ， 采 用 基于 FDR (false discovery rate) 标 
准 的 B-H 法 对 显著 性 进行 统计 校正 (Benjamini & Hochberg, 1995) 。 (5) 控制 早期 抑郁 可 
能 造成 低估 基因 和 教养 对 抑郁 的 效应 ， 本 研究 补充 考察 了 多 基因 累加 得 分 与 教养 对 T1 抑郁 
和 T2 抑郁 (不 控制 Tl 抑郁 ) 的 影响 。 采 用 SPSS 23.0 Fl R 3.3.2 进行 数据 处 理 和 统计 分 析 。 
3 结果 
3.1 三 种 基因 标记 的 基因 型 分 布 
本 研究 考察 的 COMT, DATI DRD2 基因 型 分 布 如 表 1 所 示 ， 三 个 位 点 基因 型 分 布 符 
合 Hardy-Weinberg 平衡 Cys < 1.74, df = 1, ps > 0.05) 。 三 个 基因 多 态 性 的 次 要 等 位 基因 分 
布 频率 (minor allele frequency, MAF) 均 大 于 5%。 此 外 ， 三 个 基因 之 间 不 存在 显著 关联 Os 
< 2.84, df= 4, ps > 0.58) 。 
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多 巴 胺 系统 基因 分 布 情况 


位 点 MAF H-W(X) 携带 低 活性 等 位 基因 数量 
0 1 2 
186267 79% 0.07 (p = 0.79) GG (891) GT (155) TT (6) 
1827072 23.0% 1.74 (p =0.19) TT (63) CT (357) CC (632) 
181799978 18.4% 0.01 (p = 0.93) AA (701) AG (315) GG (36) 
YE: MAF 代表 次 要 等 位 基因 分 布 频率 ，H-W 代表 Hardy-Weinberg 平衡 。 
3.2 检验 基因 标记 的 不 同 效应 
在 计算 多 基因 累加 得 分 前 ， 需 确定 各 基因 采用 线性 编码 方式 是 否 合适 、 是 否 存在 单个 基 
因 效 应 主导 多 基因 累加 效应 的 情况 ， 以 及 多 个 基因 的 效应 是 否 会 因为 编码 方式 而 相互 抵消 。 
本 研究 基于 多 巴 胺 活性 进行 显 性 和 隐 性 编码 ， 将 低 活 性 多 巴 胺 等 位 基因 编码 为 风险 / 易 感 等 
位 基因 。 参 照 既 有 研究 文献 (Stocker et al., 2017) ， 进 行 线性 基因 和 等 基因 效应 模型 检验 。 
结果 见 表 2, 分 解 模型 的 解释 率 最 高 ,不同 模 型 间 R? 改 变量 并 不 显著 。 本 研究 中 三 个 候选 基 
因 没 有 显著 偏离 线性 和 等 基因 效应 假设 。 因此， 后 续 分 析 中 ， 按 照 个 体 携带 的 低 活性 多 巴 胺 
等 位 基因 数量 进行 线性 (0、1、2) 编码 ， 并 进行 多 基因 累加 得 分 计算 (多 基因 累加 得 分 分 
fi: 1=7A, 2=43 A, 3=259 A, 4=487 A, 5=217 A, 6=39 A) 。 
R2 线性 和 等 基因 效应 模型 比较 
模型 模型 拟 合 模型 拟 合 改变 量 
R F (af) AR F (ds 
积极 教养 
模型 1: 分 解 模 型 0.373 40.15 (15, 1012) — = 
模型 2， 线 性 基因 效应 模型 0.371 66.58 (9, 1018) 0.002 0.69 (6, 1012) 
模型 3， 等 基因 效应 模型 0.365 117.63 (5, 1022) 0.006 2.13 (4, 1018) 
消极 教养 
模型 1， 分 解 模型 0.378 41.04 (15, 1012) = = 
模型 2， 线 性 基因 效应 模型 0.374 67.49 (9, 1018) 0.005 1.22 (6, 1012) 
模型 3， 等 基因 效应 模型 0.370 120.14 (5, 1022) 0.005 1.42 (4, 1018) 
母亲 教养 敏感 性 〈 补 充分 析 ) 
模型 1， 分 解 模型 0.376 40.68 (15, 1012) _ = 
模型 2， 线 性 基因 效应 模型 0.374 67.52 (9, 1018) 0.002 0.63 (6, 1012) 
模型 3， 等 基因 效应 模型 0.368 119.27 (5, 1022) 0.006 2.16 (4, 1018) 
VE: 模型 1 包括 6 个 基因 主 效应 和 6 个 基因 x 环境 交互 效应 (3 个 基因 的 线性 和 隐 性 编码 ) 以 及 控制 变量 和 教养 的 效应 ， 模 型 2 将 
同一 个 基因 的 显 性 和 隐 性 编码 回归 权重 限定 相等 (6 个 限定 ) ， 包 括 基因 主 效应 和 基因 x 环 境 交 互 效应 ;模型 3 将 3 个 基因 对 抑郁 
的 主 效应 和 基因 x 环境 交互 效应 设 定 为 等 同 (4 个 限定 ) 。 
a 模型 改变 五 值 均 不 显著 。 
3.3 各 变量 相关 分 析 结 果 


主要 研究 变量 的 平均 数 、 标 准 差 及 相关 系数 见 表 3。 多 基因 累加 得 分 与 积极 教养 、 消 极 


教养 行为 间 均 不 存在 显著 相关 ， 排 除了 
多 巴 胺 活 怕 
母亲 消极 教养 行为 显著 正 向 预测 
年 的 抑郁 症状 间 旺 上 
本 1 检验 显示 ，T2 抑郁 显著 高 于 Tl 抑郁 Caos =—5.29, p < 0.001) 


少年 间隔 一 


H 
= 


FE 越 低 ) 则 个 体 幅 是 


JOD 


抑郁 的 风险 越 大 , 母亲 积极 教养 行为 显著 负 向 预测 青 


基因 -环境 相关 的 可 能 性 


fi 


等 程度 正 相 关 , 表明 青 


少年 抑郁 。 积极 和 消极 教养 行为 i 
少年 抑郁 


。 多 基因 累加 得 分 越 高 〈 即 
少年 抑 邦 ， 
存在 显著 负 向 关联 。 
具有 一 定 的 稳定 性 , ROMP 

。 此 外 ， 与 相关 结果 相 


青 


mt 


一 致 , 独立 样本 1 检验 表明 ,母亲 积极 教养 与 抑郁 均 不 存在 显著 的 性 别 差 异 (积极 教养 : taos 


8 


= 1.46, p = 0.14; T1 抑郁 : taos) = —1.86, p = 0.06; T2 抑郁 taos) =-1.13, p = 0.26) ， 但 是 母 导 
对 男生 的 消极 教养 行为 显著 多 于 女生 (tas) =-2.54,p = 0.01) 。 
表 3 描述 统计 与 相关 系数 表 


变量 M SD 1 2 3 4 5 6 
1. 性 别 一 一 1 

2. 多 基因 累加 得 分 一 = —0.01 1 

3. 积极 教养 0.52 0.64 一 0.05 —0.01 1 

4. 消极 教养 0.30 0.43 0.08" 0.01 -0.52""* 1 

5. 教养 敏感 性 一 0.00 0.59 —0.06" —0.01 0.95"* —0.76""" 1 

6. T1 抑郁 0.22 0.24 0.06 0.07" —0.15** oir” —0.18** 1 

7. T2 抑郁 0.27 0.25 0.03 0.09** —0.12** 0.17 —0.15** 0.60°** 


TE: M 代表 了 平均 值 ，SD 代表 标准 差 ， 数 字 1~6 分 别 代表 了 变量 一 列 各 个 数字 后 面 的 变量 。 
*p<0.05; * p < 0.01; ** p < 0.001. 


34 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 抑郁 的 交互 作用 

以 T2 抑郁 为 因 变 量 ， 控 制 TL 抑郁 和 性 别 ， 以 多 基因 累加 得 分 、 母 亲 教 养 行为 以 及 基 
因 x 环 境 交 互 项 为 预测 变量 进行 分 层 回 归 分 析 。 如 表 4 所 示 ， 多 基因 累加 得 分 和 母亲 消极 教 
养 显 著 正 向 预测 青少年 抑郁 ， 但 是 母亲 积极 教养 对 抑郁 的 预测 作用 不 显著 ; 控制 Tl 抑郁 和 
性 别 后 ， 多 基因 累加 得 分 与 母亲 积极 教养 、 消 极 教养 交互 预测 青少年 抑 椰 ， 且 经 过 B-H 统 
计 校 正 后 结果 依然 显著 。 

进一步 的 简单 斜率 检验 显示 ， 在 控制 T1 抑郁 和 性 别 后 ， 在 母亲 积极 教养 行为 较 低 时 ， 
具有 较 高 多 基因 累加 得 分 的 青少年 比 具 有 较 低 多 基因 累加 得 分 的 青少年 抑郁 水 平 更 高 (pb = 
0.03, t= 2.89, p= 0.004) ， 而 在 母亲 积极 教养 行为 较 高 时 ， 不 同 多 基因 累加 得 分 的 青少年 抑 
郁 水 平 无 显著 差异 (b = -0.002, t=-0.21, p = 0.84) 。 与 此 类 似 ， 在 母亲 消极 教养 行为 较 高 
时 ， 有 具有 较 高 多 基因 累加 得 分 的 青少年 比 具 有 较 低 多 基因 累加 得 分 的 青少年 抑郁 水 平 更 高 
Cb = 0.03, t= 3.18, p = 0.002) ， 而 在 母亲 消极 教养 行为 较 低 时 ， 不 同 多 基因 累加 得 分 的 
少年 抑郁 水 平 无 显著 差异 Cb = -0.003,t= -0.42, p=0.68) ( 详 见 图 1) 。 


a 


LLL 
iu 


4 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 抑郁 的 影响 


AR? b SE B p 
模型 1 
性 别 0.359 一 0.00 0.01 —0.01 0.76 
早期 抑郁 0.64 0.03 0.59 <0.001 
多 基因 累加 得 分 (PS) 0.003 0.01 0.01 0.05" 0.04 
积极 教养 一 0.01 0.01 一 0.03 0.24 
PSx 积 极 教养 0.003 一 0.01 0.01 —0.05* 0.03 
模型 2 
性 别 0.359 —0.02 0.01 —0.01 0.66 
早期 抑郁 0.63 0.03 0.59*** <0.001 
多 基因 累加 得 分 (PS ) 0.01 0.01 0.01 0.05" 0.05 
消极 教养 0.02 0.01 0.07"* 0.01 
PSx 消 极 教养 0.004 0.01 0.01 0.06" 0.01 
模型 3 〈 补 充分 析 ) 
性 别 0.359 —0.00 0.01 —0.01 0.71 
早期 抑郁 0.63 0.03 0.59 <0.001 
多 基因 累加 得 分 (PS) 0.01 0.01 0.01 0.05" 0.04 
教养 敏感 性 一 0.01 0.01 ~0.05 0.048 
PSx 教 养 敏 感性 0.01 0.02 0.01 0.07™ 0.006 


注 : *p<0.05; p<0.01;" p < 0.001. 


oe 高 多 基因 累加 得 分 高 多 基因 累加 得 分 03 。 一 一 高 基因 累加 得 分 
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图 1 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 抑郁 的 影响 

如 表 5 所 示 , 本 研究 分 别 对 积极 教养 和 消极 教养 模型 进行 再 参数 化 回归 , 检验 上 述 交 互 
作用 符合 何 种 理论 模型 (“ 素 质 - 压 力 ” 模 型 与 “不 同 易 感 性 ”模型 ) 。 结 果 显示 ， 无 论 是 积 
极 教养 还 是 消极 教养 模型 ， 多 基因 累加 得 分 x 教 养 的 交互 作用 均 符 合 “ 素 质 -压力 ”模型 ， 详 
细 统 计 指标 分 析 如 下 : G) “素质 -压力 ”模型 本 身 拟 合 良好 (R?=0.37,p<0.001) ; (2) 
在 “不 同 易 感性 ”模型 中 ， 交 叉 点 C 的 点 估计 值 及 其 95% 置 信 区 间 超 出 了 教养 行为 的 取 值 
范围 (积极 教养 -4.22~ 1.52; 消极 教养 : -1.31~ 3.76; 教养 敏感 性 : 4.16~1.60) ， 这 表 
明基 因 x 环 境 交 互 模式 不 符合 “不 同 易 感性 ”假设 。 G) 尸检 验 显 示 尽 管 “ 不 同 易 感性 ” 模 
型 比 “ 素 质 -压力 ”模型 增加 一 个 参数 (C) 估计 ， 但 “不 同 易 感 性 ”模型 的 解释 率 并 没有 显 
著 高 于 “素质 -压力 ”模型 (AR2s= 0.00,ps> 0.05) 。 (4) “素质 -压力 ”模型 的 AIC 和 BIC 
指标 均 小 于 “不 同 易 感性 ”模型 。 综 上 ， 各 参数 均 表明 , “素质 -压力 ”模型 对 数据 的 拟 合 优 
于 “不 同 易 感性 ”模型 ， 因 此 本 研究 中 多 巴 胺 系统 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 
抑郁 的 交互 作用 更 加 符合 “素质 -压力 ”模型 。 

表 5 多 基因 累加 分 与 父母 教养 行为 对 青少年 抑郁 交互 作用 的 再 参数 化 回归 模型 检验 


(VY 积极 教养 行为 消极 教养 行为 教养 敏感 性 总 分 (补充 分 析 ) 
= “不 同 易 感性 ” “素质 压力 ” “不 同 易 感性 ” “素质 -压力 ” “不 同 易 感性 ” “素质 -压力 ” 
K. Bı —0.00(0.01) —0.00(0.01) —0.01 (0.01) —0.01 (0.01) —0.00(0.01) —0.00(0.01) 

© By 0.64(0.03)"** 0.63(0.03)"** 0.63(0.03)"** 0.62(0.03)** 0.63(0.03)"** 0.63(0.03)"** 

Bs —0.01(0.01) —0.01(0.01) 0.02(0.01)* 0.02(0.01)* ~0.01(0.01)* ~0.01(0.01)* 
By —0.01(0.01)* —0.01(0.003)* 0.01(0.01)* 0.01(0.003)* ~0.02(0.01)* —0.01(0.003)* 
C 0.88(0.58) 1.52(—) 一 0.77 —1.31(—)' 0.69 (0.41) 1.60(—)" 
95% CI of C [0.03, 12.31] = [-3.79, 0.01] — [0.03, 2.69] — 
R 0.37 0.37 0.37 0.37 0.37 0.37 
F(df) 117.63(5, 1022)"**  146.93(4, 1023)"** 117.97(5, 1022)*** 150.04(4, 1023)""* 119.25(5, 1022)"* 148.48(4,1023)™ 
F a vs. b (df) = 0.64(1, 1022) 至 0.71(1, 1022) = 1.92(1, 1022) 
AIC 一 554.82 一 556.18 一 562.81 一 564.09 -560.06 -560.13 
BIC 一 520.27 一 526.57 一 528.26 一 534.48 -525.52 -530.52 


YE: 模型 了 am = Bot Byx PES + B2x 早 期 抑郁 + Bs x (Xas — C)+ Bax (Xar 一 C)x 多 基因 累加 得 分 +E; CI 代表 置信 区 间 ; F 
a vs. b 代表 “不 同 易 感性 ”与 “素质 -压力 ”模型 之 间 的 下 检验 ， 以 分 析 减 少 一 个 估计 参数 时 模型 的 AR? 是 否 存在 显著 变化 ，a 
代表 该 参数 限定 在 指定 数值 上 ， 如 C = 1.52、 一 1.31、1.60 表示 将 交叉 点 C 固定 在 积极 教养 的 最 大 值 1.52、 消 极 教养 的 最 小 值 一 
1.31 处 、 母 亲 教 养 敏感 性 〈 积 极 和 消极 教养 合并 分 ) 的 最 大 值 1.60 处 。 

* p<0.05; * p < 0.01; ** p< 0.001. 
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3.5 敏感 性 分 析 


为 验证 上 述 研 究 结果 的 可 靠 性 ， 本 研究 进行 了 敏感 性 分 析 。 首 先 ， 对 任何 单 基因 与 母亲 
积极 教养 、 消 极 教养 行为 对 抑郁 的 交互 作用 进行 回归 分 析 。 结 果 显 示 在 B-H 校正 后 , BR DATI 


基因 与 积极 教养 行为 的 交互 作用 显著 外 (b = 0.05, t =-2.78, p= 0.006) ， 任 何 单 基因 与 教 
养 的 交互 作用 均 不 显著 (|b|s < 0.04, ltls = -2.27, ps > 0.05) 。 这 也 表明 多 基因 累加 得 分 研究 


大 
= 
ai 
at 
At 
发 
Eo 
w 
>H 
a 


著 性 方面 比 单 基 


因 研 究 更 具 优 势 。 


基因 累加 得 分 作为 连续 变量 进行 分 析 ， 由 于 携带 1、2 和 6 个 低 多 


果 ， 本 研究 对 多 基因 累加 得 分 重新 编码 进行 
青少年 编码 为 低 分 组 , 将 多 基 


因 的 人 数 较 少 , 可 能 影响 上 述 交 互 作用 结果 的 可 靠 性 。 为 进一步 验证 上 述 


分 为 5、6 的 青少年 编码 为 高 分 组 。 


因 累 加 得 分 为 4 的 青少年 编码 为 中 等 得 分 组 , 将 多 基因 累加 得 
分 组 回归 分 析 结果 与 前 述 结果 相 一 致 ， 相 比 低 分 组 和 中 


分 组 回归 检验 。 将 多 基因 累加 得 分 为 1、2、3 的 


等 得 分 组 ， 高 分 组 的 青少年 对 母亲 积极 教养 (b= -0.04, t= -2.33,p=0.02〉 和 消极 教养 更 加 
敏感 (b= 0.04, t=2.65,p=0.01) ( 详 见 图 2) 。 
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图 2 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 抑郁 影响 的 敏感 性 分 析 


第 三 ， 为 了 排除 基因 x 基因 交互 效应 的 影响 ， 本 研究 对 G x GxGxE 效 应 进行 了 分 析 。 


需要 指出 的 是 ， 采 用 线性 基 


因 编 码 (0、 


制 无 法 满足 统计 分 析 条 件 。 本 研究 痊 代 性 


0=T, 1 = CC; DRD2:0=AA,1=G) i 


1, 2) 的 方式 会 导致 每 个 基因 型 组 合 单元 格 人 数 限 
的 采用 了 二 分 编码 方式 (COMT7: 0=GG, 1=TD477: 


行 分 析 ， 结 果 显 示 任 何 GxG，GxGxG，GxGx 


E，GxGxGxE 交 互 效应 均 不 显著 〈lbls< 0.22, ltls < 0.90, p>0.37) ， 在 一 定 程度 上 排除 
了 潜在 基因 x 基因 效应 的 干扰 。 


第 四 ， 参 照 既 有 研究 文献 (Cao et al., 2018) ， 将 消极 教养 和 积极 教养 合并 后 的 母亲 教 


养 敏 感性 总 分 作为 环境 指标 , 重复 上 述 统计 分 析 , 结果 与 单独 进行 积极 和 消极 教养 行为 的 模 


式 一 致 ， 只 是 交互 项 的 显著 
第 五 ， 补 充 考 察 了 多 基 
向 效应 (不 控制 Tl 抑郁 ) 。 


结果 显示 ， 


性 和 效应 量 更 大 参见 表 2、 表 4 和 表 5 的 补充 分 析 部 分 ) 。 
因 累 加 得 分 和 教养 对 Tl 抑郁 的 同时 效应 ， 以 及 对 T2 抑郁 的 纵 


多 基因 与 教养 的 主 效应 和 显著 性 提高 。 但 是 ， 除 基 


因 与 消极 教养 对 T2 抑郁 的 交互 作用 显著 外 ， 其 他 基因 x 环境 均 不 显著 〈 见 表 6 和 表 7) 。 


表 6 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 T1 抑郁 的 影响 


AR? b SE Bp p 

模型 1 

性 别 0.003 0.02 0.01 0.05 0.09 

多 基因 累加 得 分 (PS) 0.03 0.01 0.01 0.07" 0.03 

积极 教养 -0.03 0.01 —0.15*** < 0.001 

PSx 积 极 教养 0.001 0.01 0.01 0.03 0.27 
模型 2 

性 别 0.003 0.02 0.01 0.05 0.14 

多 基因 累加 得 分 (PS ) 0.03 0.01 0.01 0.07 0.02 

消极 教养 0.04 0.01 0.17” < 0.001 

PSx 消 极 教养 < 0.001 0.00 0.01 -0.01 0.87 
模型 3 

性 别 0.003 0.02 0.01 0.05 0.13 

多 基因 累加 得 分 (PS ) 0.04 0.01 0.01 0.07" 0.03 

教养 敏感 性 —0.04 0.01 -0.17"" < 0.001 

PSx 教 养 敏感 性 0.001 0.01 0.01 0.03 0.35 


yÈ: “p<0.05; p<0.01 


;™" p< 0.001. 


表 7 多 基因 累加 得 分 与 母亲 教养 行为 对 青少年 T2 抑郁 的 影响 〈 未 控制 早期 抑郁 ) 


AR? b SE B p 
模型 1 
性 别 0.001 0.01 0.01 0.03 0.37 
多 基因 累加 得 分 (PS) 0.02 0.02 0.01 0.09** 0.002 
积极 教养 —0.03 0.01 -0.12""* < 0.001 
PSx 积 极 教养 0.001 -0.01 0.01 ~0.03 0.30 
模型 2 
性 别 0.001 0.01 0.01 0.02 0.51 
多 基因 累加 得 分 (PS) 0.04 0.02 0.01 0.09™ 0.002 
消极 教养 0.04 0.01 0.17°* < 0.001 
PSx 消 极 教养 0.003 0.01 0.01 0.06" 0.05 
模型 3 
性 别 0.003 0.01 0.01 0.02 0.45 
多 基因 累加 得 分 (PS ) 0.03 0.02 0.01 0.09** 0.002 
教养 敏感 性 —0.04 0.01 -0.15*"* < 0.001 
PSx 教 养 敏感 性 0.002 -0.01 0.01 —0.05 0.10 
YE: *p<0.05;* p< 0.01; ** p < 0.001. 
4 讨论 


本 研究 采用 多 基因 累加 得 分 研究 范式 ， 以 三 种 多 巴 胺 系统 基因 (COMT 基因 、D477 基 


因 和 DRD2 基因 


) 为 遗传 指标 ， 考 察 了 母亲 教养 行为 与 多 基因 累加 得 分 对 青少年 抑 邦 


的 影 


啊 。 结 果 显 示 控 制 早期 抑郁 后 ， 多 基因 累加 得 分 显著 正 向 预测 青少年 抑郁 ; 多 巴 胺 系统 


因 累 加 得 分 与 母 杀 积极 、 消 极 教养 均 交 互 预测 青少年 抑郁 ,其 交互 模式 均 符 合 “ 素 质 - 压 


模型 。 具 体 表现 为 ， 在 低 积 极 或 高 消极 教养 环境 中 ， 多 基因 累加 得 分 较 高 的 青少年 比 多 
累加 得 分 较 低 的 青少年 抑郁 水 平 更 高 , 但 是 在 高 积极 或 低 消极 教养 环境 中 , 携带 不 同 多 
累加 得 分 的 青少年 抑郁 水 平 不 存在 显著 差异 。 

与 以 往 研究 (e.g., Wang etal.,2016) 相 一 致 ， 本 研究 发 现 母 亲 消 极 教养 行为 是 导致 
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年 抑郁 的 风险 因素 。 但是， 母亲 积极 教养 对 青少年 抑郁 的 影响 未 达到 显著 性 水 平 。 这 可 外 
由 于 母亲 积极 教养 和 消极 教养 行为 对 青少年 抑郁 的 影响 存在 差异 。 一 方面 , 积极 和 消极 
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行为 对 青少年 抑郁 的 作用 程度 存在 差异 ， 如 Dallaire 等 (2006) 的 研究 发 现 ， 消 极 教养 行为 
对 儿童 青少年 抑郁 的 预测 力 更 强 。 另 一 方面 , Smokowski，Bacallao，Cotter 和 Evans (2015) 
研究 进一步 揭示 , 不 同类 型 的 教养 行为 对 青少年 适应 结果 的 预测 作用 存在 领域 特殊 性 , 积极 
教养 可 能 与 积极 心理 特征 的 关联 更 为 密切 , 而 对 消极 适应 结果 的 影响 较 小 , 而 消极 教养 行为 
则 是 焦虑 、 抑 郁 等 消极 适应 结果 的 关键 风险 因素 。 更 重要 的 是 , 积极 和 消极 教养 对 青少年 抑 
郁 的 作用 路 径 可 能 也 存在 差异 ， 如 Luebbe 和 Bell (2014〉 的 研究 发 现 消极 家 庭 氛围 (包括 
消极 教养 行为 ) 和 积极 家 庭 氛围 《包括 积极 教养 行为 ) 分 别 通过 消极 情绪 性 和 积极 情绪 性 间 
接 影响 青少年 抑郁 。 由 此 , 未 来 研究 可 以 增加 积极 适应 结果 考察 教养 与 适应 结果 间 的 领域 特 
殊 性 问题 。 

本 研究 发 现 个 体 携带 的 多 基因 累加 得 分 越 高 (多 巴 胺 系统 功能 越 低 ) 则 其 抑郁 水 平 越 高 。 
这 与 单 胺 缺陷 假说 相 一 致 , 多 巴 胺 系统 功能 的 缺陷 会 增加 个 体 纵 患 抑郁 的 风险 .更 重要 的 是 ， 
本 研究 为 抑郁 的 多 基因 遗传 结构 提供 了 支持 。 但是, 本 研究 中 多 基因 联合 效应 的 解释 率 仍然 
为 1% 左右 。 这 可 能 提示 了 抑郁 具有 复杂 的 产生 和 发 展 机 制 ， 本 研究 所 选取 的 多 巴 胺 系统 多 
基因 累加 得 分 也 仅 能 反映 其 致 病因 素 的 冰山 一 角 , 研究 者 应 该 更 加 慎重 地 看 待 这 一 结果 。 此 
Sh, 研究 者 应 该 深入 探索 多 基因 影响 抑郁 的 作用 机 制 , 如 “多 基因 累加 得 分 x 环境 ”研究 等 。 

从 交互 作用 上 来 看 , 在 控制 早期 抑郁 后 , 母亲 教养 行为 与 多 巴 胺 系统 多 基因 累加 得 分 交 
互 预测 青少年 后 继 抑郁 , 并 且 进 一 步 的 分 析 发 现 ， 只 有 在 消极 环境 下 (包括 低 积 极 教养 和 高 
消极 教养 ) ,多 基因 累加 得 分 较 高 的 青少年 比 得 分 较 低 的 个 体 具有 更 高 的 抑郁 风险 。 再 参数 
化 回归 分 析 进 一 步 显示 这 一 交互 作用 模式 更 符合 “素质 -压力 ”模型 。 与 此 结果 相 一 致 ， 前 
述 关 于 腹 侧 纹 状 体 活性 的 研究 显示 , 个 体 携带 的 低 活 性 多 巴 胺 系统 等 位 基因 越 多 , 则 在 面 对 
外 界 刺激 时 表现 出 更 低 的 纹 状 体 反应 性 (Nikolova et al., 2011) 。 有 关 奖 赏 和 惩罚 敏感 性 的 
研究 发 现 , 对 奖赏 刺激 的 敏感 性 与 纹 状 体 激活 水 平 呈 正 向 关联 ,而 对 惩罚 的 敏感 性 与 纹 状 体 
激活 程度 呈 负 向 关联 (Kim, Yoon, Kim, & Hamann, 2015) 。 这 些 研究 结果 提示 ， 携 带 低 活性 
多 巴 胺 等 位 基因 数量 越 多 的 个 体 具 有 较 低 的 纹 状 体 激活 水 平 , 对 社会 奖赏 刺激 (如 积极 的 教 
养 环境 ) 的 敏感 性 较 低 从 而 无 法 在 积极 的 环境 中 获 益 , 但 是 其 对 消极 环境 刺激 的 敏感 性 较 高 
从 而 增加 其 抑郁 的 风险 。 需要 指出 的 是 , 这 一 交互 作用 模式 可 能 也 仅仅 描绘 了 青少年 早期 阶 
段 多 基因 x 环 境 的 交互 作用 ， 当 前 研究 者 更 倾向 于 采用 动态 的 视角 看 待 “ 素 质 - 压 力 ” 和 “不 
同 易 感 性 ”模型 ( 王 美 萍 ,2015) 。 壁 如, 研究 发 现 多 基因 累加 得 分 CBDNF 基因 和 5-HTTLPR 
基因 ) 与 家 庭 环 境 质 量 对 抑郁 的 交互 作用 模式 在 青少年 早期 (15 岁 前 ) 符合 “不 同 易 感 性 ” 
假说 ， 而 在 青少年 中 晚期 则 符合 “素质 -压力 ”模型 (Dalton, Hammen, Najman, & Brennan, 
2014) 。 由 此 ， 未 来 研究 应 该 补充 其 他 年 龄 阶段 的 被 试 或 采用 追踪 设计 ， 深 入 探索 多 基因 
环境 交互 作用 模式 是 否 具 有 发 展 动态 性 。 

虽然 本 研究 中 多 基因 累加 得 分 及 其 与 环境 的 交互 作用 解释 率 仍 然 很 低 〈1%) ， 但 敏感 
性 分 析 显 示 , 与 单 基因 相 比 , 多 基因 累加 得 分 在 提高 遗传 解释 率 和 显著 性 方面 仍然 具有 一 定 
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的 优势 。 值得 指出 的 是 ,这 并 不 意味 着 研究 者 可 以 为 扩大 遗传 效应 而 随意 增加 基因 位 点 数量 
来 计算 多 基因 累加 得 分 。 事 实 上 ， 即 便 累 加 了 众多 SNP 基因 位 点 ， 也 无 法 解释 抑郁 全 部 的 
遗传 机 制 ， 甚 至 对 遗传 率 的 提升 甚 微 。 如 Peyrot “ (2014) 采用 全 基因 组 关联 扫描 (GWAS) 
方法 ， 分 别 将 150~32870 个 SNP 位 点 纳入 多 基因 累加 得 分 计算 ， 其 多 基因 累加 得 分 及 其 与 
童年 期 虐待 的 交互 作用 解释 率 仅 为 0.09%~0.90%。 但 是 , 较 小 的 遗传 解释 率 并 不 意味 着 它 是 
可 以 忽视 的 (Evans, 1985) 。 因 此 ,研究 者 不 能 盲目 追求 解释 率 的 提升 而 滥用 遗传 指标 。 由 
此 ,如 何 选 择 多 基因 累加 指标 成 为 研究 者 关注 的 重点 , 伴随 着 越 来 越 多 新 功能 性 基因 标记 的 
发 现 ， 这 一 问题 变 得 尤为 突出 。 对 此 ，Vrshek-Schallhorn “ (2015) 指出 为 了 避免 数据 驱 
和 小 用 的 风险 , 在 进行 多 基因 累加 效应 研究 时 最 好 选择 同一 神经 递 质 系 统 的 基因 而 不 是 跨 系 
统 的 基因 , 这 样 有 助 于 探索 多 基因 与 内 表 型 之 间 的 关联 , 进而 确定 抑郁 发 生 的 内 部 作用 机 制 。 
此 外 ， 多 基因 累加 得 分 解释 率 的 限制 使 得 研究 者 对 探索 “缺失 遗传 率 ” 表 现 出 极 大 的 兴 
E 趣 。Manolio 等 〈2009) 指出 遗传 率 缺 失 的 原因 可 能 包括 : 第 一 ， 大 量具 有 极 小 效应 的 基因 
变 体 尚未 被 发 现 。 因 此 ,采用 GWAS 和 GWEIS( genome-wide by environment interaction studies ) 
研究 方法 探索 潜在 的 功能 性 基因 及 潜在 的 基因 x 环境 交互 作用 仍然 是 未 来 研究 的 重点 方向 。 
第 二 ， 研 究 者 应 该 跳出 “常见 疾病 ， 常 见 变 体 〈common disease, common variant) ”的 固有 
假设 ， 探 索 遗 传 效 应 较 大 的 罕见 基因 变 体 〈 在 人 群 中 的 比例 低 于 5$% 的 基因 变 体 ) 。 第 三 ， 
现 有 研究 很 少 测量 DNA 表达 机 制 ， 限 制 了 其 对 表 型 的 遗传 解释 率 。 甲 基 化 、 组 蛋白 修饰 等 
表 观 遗传 机 制 可 以 在 不 改变 DNA 结构 的 情况 下 影响 基因 - 表 型 间 的 关联 。 因 此 , 考察 多 基因 
累加 得 分 与 基因 表达 水 平 的 变化 , 可 用 于 探寻 其 遗传 解释 率 的 缺失 并 更 进一步 剖析 遗传 因素 
的 潜在 生物 机 制 。 第 四 ， 多 基因 累加 得 分 未 能 检测 基因 x 基因 交互 作用 。 在 本 研究 中 ， 受 单 
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元 格 人 数 限 制 ， 采 用 二 分 编码 方式 对 基因 x 基 因 交 互 效应 进行 了 检验 并 在 一 定 程度 上 证 实 了 
= 多 基因 累加 的 可 靠 性 。 但 是 , 在 未 来 多 基因 累加 研究 中 尚 需 对 不 同 基因 间 的 交互 效应 保持 说 
慎 。 正 如 Plomin，DeFries，McClearn 和 McGuffin (2001) 指出 多 基因 遗传 模型 除了 需要 考 
虑 不 同等 位 基因 效应 间 的 累加 效应 外 , 还 需要 考虑 不 同 基因 座 的 效应 彼此 间 有 相互 作用 的 可 
能 性 。 

本 研究 采用 纵向 设计 控制 早期 抑郁 , 考察 多 基因 累加 得 分 与 母亲 教养 对 青少年 抑 邦 的 影 
响 。 采 用 这 一 设计 的 主要 原因 是 青少年 抑郁 与 母亲 教养 间 通 常 存在 双向 关系 〈Hamza & 
Willoughby, 2011) , 即 母亲 教养 不 仅 影 响 青少年 抑郁 , 而 且 也 会 受到 青少年 抑郁 情绪 的 影响 。 
由 此 ,如果 不 控制 青少年 早期 抑郁 , 教养 行为 对 青少年 后 继 抑郁 的 影响 本 质 上 可 能 也 混杂 了 
早期 抑郁 对 后 继 抑 郁 的 影响 。 因 此, 控制 早期 抑郁 有 助 于 降低 母亲 教养 行为 是 抑郁 症状 的 结 
果 而 不 是 原因 的 可 能 性 〈Chen, Li, & McGue, 2013) 。 虽 然 控 制 早 期 抑郁 会 带 来 一 些 优 势 ， 
但 也 可 能 会 造成 基因 和 环境 的 效应 量 和 显著 性 降低 。 此 外 , 综合 考虑 控制 早期 抑郁 与 不 控制 
早期 抑郁 的 结果 , 在 理解 本 研究 结果 时 需要 注意 的 是 , 多 基因 累加 得 分 与 母亲 教养 的 交互 作 
用 不 能 预测 Tl 和 T2 时 间 点 上 绝对 水 平 的 抑郁 症状 ,本 质 上 是 预测 了 T1 到 T2 抑郁 的 变化 。 


14 


— 


在 Petersen (2012) 的 研究 中 也 发 现 了 类 似 的 结果 模式 ， 即 G x 下 不 能 预测 基线 水 平 的 焦 
虑 抑郁 症状 ,但 能 够 显著 预测 焦虑 抑郁 的 发 展 变化 , 尤其 是 能 够 解释 抑郁 随 年 龄 的 变化 趋势 。 
由 此 ， 即 使 未 能 发 现 遗 传 与 环境 交互 效应 对 抑郁 水 平 的 显著 预测 作用 ， 也 不 意味 着 GxE 与 
抑郁 发 展 无 关 。 未 来 研究 应 该 更 深入 的 考察 遗传 与 环境 对 抑郁 及 其 发 展 变化 的 影响 以 更 深入 
地 揭示 此 问题 。 

需要 注意 的 是 , 本 研究 也 具有 一 些 局 限 性 ,第 一 , 本 研究 缺乏 重复 验证 性 样本 。 近年 来 ， 
伴随 着 基因 x 环境 交互 作用 研究 分 歧 的 增多 ， 越 来 越 多 的 研究 者 对 遗传 研究 的 可 重复 性 问题 
提出 了 质疑 ， 并 且 指出 应 该 采用 外 部 验证 样本 进行 重复 检验 (Christ, Schwartz, Stoltenberg, 
Brauer, & Savolainen, 2018; Dick et al., 2015) 。 但 是 本 研究 未 能 具有 相同 测量 工具 及 相似 被 
试 特征 的 外 部 验证 样本 ， 其 结果 的 可 重复 性 有 竺 未 来 研究 进行 验证 。 第 二 ， 多 巴 胺 系统 基因 
众多 , 本 研究 仅 采用 了 三 个 基因 多 态 性 考察 了 多 基因 遗传 基础 ,未 来 研究 尚 需 要 考察 其 他 候 
选 基因 多 态 性 以 丰富 该 领域 的 研究 。 第 三 ， 本 研究 中 消极 教养 的 信 度 指标 较 低 ， 这 可 能 与 消 
极 教养 的 测量 题目 (4 个 ) 较 少 有 关 。 虽 然 本 研究 的 验证 性 分 析 也 保证 了 其 结构 效 度 ， 但 仍 
需要 谨慎 的 看 待 本 研究 结果 。 此 外 ， 本 研究 仅 考 察 了 母亲 积极 和 消极 教养 行为 ,， 未 对 父亲 教 
养 行 为 进行 测量 。 长 期 以 来 ， 在 分 子 遗 传 研究 领域 , 研究 者 往往 忽视 了 父亲 、 母 亲 教 养 行 为 
的 差异 性 及 相互 影响 。 据 我 们 所 知 ， 目 前 仅 有 一 项 研究 比较 了 多 基因 累加 得 分 与 父亲 、 母 亲 
教养 行为 交互 作用 的 不 同 但 未 发 现 显著 差异 〈Stockeretal.,2017) 。 但 是 ， 也 有 
亲 和 母 亲 教 养 行为 交互 影响 青少年 抑郁 (Quach, Epstein, Riley, Falconier, & Fang,2015) 。 由 
此 ， 未 来 研究 应 该 考察 父亲 、 母 亲 教 养 行为 与 遗传 因素 交互 作用 的 差异 ， 并 且 关注 父亲 和 和 母 
亲 教 养 行 为 间 的 相互 影响 。 最 后 ， 本 研究 被 试 来 自 城市 常态 群体 青少年 ， 该 研究 结果 能 否 适 
于 临床 样本 或 重 性 抑郁 以 及 社会 经 济 地 位 较 低 的 青少年 群体 尚 需 未 来 研究 的 验证 
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The Influence of Dopaminergic Genetic Variants and Maternal 
Parenting on Adolescent Depressive Symptoms: A Multilocus 
Genetic Study 
CAO Yanmiao; ZHANG Wenxin 
(School of Psychology, Shandong Normal University, Jinan 250014, China) 
Abstract 

For decades, there is increasing evidence for the importance of single-gene by environment 
interactions (G x E) in understanding the etiology of depression. However, several concerns have 
been raised about the ignoring the polygenic traits of depression when conducting G x E research 
using single loci. Within this context, the multilocus genetic profile score (MGPS) have recently 
emerged as an approach of capturing polygenic nature across multiple genes. In line with the 
monoamine deficiency hypothesis, recent research has begun to show that the combined effects of 
multiple dopaminergic genetic variants are stronger than the influence of any single gene examined 
in isolation. Additionally, genes related to the functioning of the dopaminergic system, which 
coordinates individual’s response to stress. However, existing G x E research has largely focused 
on adverse family environments (i.e., maltreatment, maternal unresponsiveness) and to a lesser 
extent on positive environment, such as positive parenting. Therefore, the present study aimed to 
examine the interaction between dopaminergic genetic variants and maternal parenting on 


adolescent depressive symptoms, by adopting the approach of multilocus genetic profile score. 


Participants were 1052 mother—offspring (adolescents mean age 12.31 + 0.37 years old at the 
first time point, 50.2% females) dyads recruited from the community. Youth completed assessments 
twice with an interval of one year. Saliva samples, self-reported depressive symptoms and mother- 
reported parenting were collected. All measures showed good reliability. Genotyping in three 
dopaminergic genes were performed for each participant in real time with MassARRAY RT software 
version 3.0.0.4 and analyzed using the MassARRAY Typer software version 3.4 (Sequenom). To 
examine whether multilocus genetic profile score moderates the effects of parenting on adolescent 
depressive symptoms and whether this potential moderating effect act in a diathesis—stress or 
differential susceptibility manner, hierarchical regression analyses were conducted. We also tested 
above questions by recoding into categorical variables and re-conducted analyses. 

The results found that multilocus genetic profile score was a significant risk factor of depression, 
with higher dopamine genetic risk scores (indicating lower dopaminergic neurotransmission) 
predicting higher levels of depression. After controlling for gender and prior depressive symptoms, 
the G x E effect with positive and negative parenting were also significant, suggesting that G x E 
interaction significantly predicted change in depression level between Time 1 and Time 2. 
Specifically, adolescents with higher MGPS exhibited higher risk for depression when encountered 
with lower levels of positive parenting and higher levels of negative parenting, compared to their 
counterparts with lower MGPS. The results support the diathesis-stress model and highlight the 
complex ways that genes and environment interact to influence development. 

These finding underscores complex polygenic underpinnings of depression and lends support 
for the mulitlocus genetic profile scores—environment interactions implicated in the etiology of 
depressive symptoms. 

Key words adolescent depression, maternal parenting, dopamine, multilocus genetic profile score, 


gene-environment interaction 


